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Abstract 

We perform a comparative study of the ratio of lattice (Wilson fcrmion) renormalization 
constants Zp/Zs, obtained non-perturbatively from the RI/MOM renormalization conditions 
and from Ward Identities of on- and off-shell Green's functions. The off-shell Ward Identity 
used in this work relies on correlation functions with non-degenerate quark masses. We find 
that, due to discretization effects, there is a 10%-15% discrepancy between the two Ward Iden- 
tity determinations at current bare couplings (j3 values). The RI/MOM result is in the same 
range and has a similar systematic error of 10-15%. Thus, contrary to a previous claim, the 
contamination of the RI/MOM result from the presence of a Goldstone pole at renormalization 
scales /i ~ a -1 is subdominant, compared to finite cutoff effects. 



1 Introduction 



The Non-Perturbative (NP) renormalization technique proposed in Q has been successfully applied 
to compute renormalization constants of two-fermion ||, || and four-fermion Q operators with 
Wilson fermions. The standard implementation consists in imposing the RI/MOM renormalization 
conditions on conveniently defined amputated projected Green functions computed between off- 
shell quark states of operators of interest, evaluated numerically at fixed momentum and lattice 
spacing (fixed lattice coupling j3). In order to control discretization effects, one should work at 
renormalization scales fi well below the UV cutoff (i.e. fi <C C(a -1 )). Moreover, any contribution 
in the renormalization constants due to spontaneous chiral symmetry breaking terms in the Green's 
functions has to be avoided to ensure that the non-perturbatively renormalized operators satisfy 
the standard QCD axial Ward Identities (abbreviated as Wis), i.e. Aqcd <S M- Thus, the existence 
of a "renormalization window" Kqcd ^ i 1 *C @( al ) is essential to the reliability of the RI/MOM 
scheme in practical simulations, as was pointed out in ref. pj]. We note in passing that, as discussed 
in refs. || the lower bound can be relaxed by matching renormalization conditions between a 
coarse and a fine lattice in the spirit of ref. || ) . 

This has not always been the case at current values of the lattice coupling. In off-shell cor- 
relation functions, at high enough momenta, the dominant contributions are perturbative. For 
example, the leading non perturbative contribution of the pseudoscalar correlation function can be 
easily obtained from the Operator Product Expansion (OPE) of the quark propagator, which has 
been worked out in refs. Q. Using either the lattice || or the Sum Rules j|] determination of the 
chiral condensate and quark mass values typical of present-day simulations, we estimate that the 
non perturbative contribution to the pseudoscalar correlation function at fi ~ 2 GeV is at most 2%. 
Therefore, the existence of a "renormalization window' has been checked in |2|, [k| by comparing 
the p 2 dependence of the Green's functions to the logarithmic one predicted by continuum pertur- 
bation theory []. In this sense, at (3 = 6.0, 6.2 and 6.4, the scalar density renormalization constant 
Zs is well behaved in a large ^-range, whereas there is no clear evidence of such a satisfactory 
renormalization window for the RI/MOM renormalization constant of the pseudoscalar density Zp 
(see ref. [0 for details). 

Recently an alternative interpretation of the observed discrepancy between the RI/MOM and 
perturbative Zp results has been proposed in ref. 11]. It is suggested that the p 2 dependence of the 
pseudoscalar correlation function is due to large contributions of the spontaneous chiral symmetry 
breaking terms, which persist at scales used in present day simulations (i.e. fi ~ a -1 ). In other 
words, the lower bound of the renormalization window is not adequately satisfied. In order to 
test this proposal, the authors of ref. [11] analyze the pseudoscalar correlation function, taking 
into account the leading OPE non-perturbative contribution and ignoring discretization effects. 
They find that "a stricking and unexpected feature of the lattice data is the very large size of the 
Goldstone boson contribution to the pseudoscalar vertex" . 

In the present work we critically examine the claim of ref. Jl l| by comparing the RI/MOM 



1 Note that continuum perturbation theory converges much better than the lattice one. Moreover, higher orders 
have been calculated in the former case, whereas only 1-loop results are available in lattice perturbation theory. 
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determination of the ratio Zp/Zg to those obtained from two Wis at /3 = 6.2,6.4 (the quenched 
approximation and Wilson fermions are implied throughout). WI results are unaffected by Gold- 
stone pole contamination but are subject to discretization effects. Thus any discrepancy between 
two WI results gives an estimate of 0{a) effects, whereas an eventually big discrepancy between 
the WI results and the RI/MOM ones would signal a Goldstone pole contamination in Zp. We 
find that individual WI results for ZpjZg are characterized by errors of about 5% but they are 
discrepant to each other at the level of 10-15% due to 0(a) effects. The RI/MOM result obtained 
at [i ~ 2 GeV is compatible with the WI ones but has a larger error of about 10-15%. Thus, 



contrary to the claim of ref. 11], we conclude that the RI/MOM determination of ZpjZg (and 
subsequently of Zp) is affected by large discretization errors while any contamination from Gold- 
stone pole contributions cannot be discerned. Analogous conclusions, based on an analysis of the 
form factors of the improved quark propagator, have been drawn in ref. |l<q| . 

Finally we wish to stress that the WI proposed in this paper (based on off-shell correlation 
functions and non-degenerate quark masses) has given a very stable estimate of ZpjZ$- In view of 
this, we prefer to evaluate Zp from the WI result of ZpjZg and the RI/MOM value of Z$- 



2 The RI/MOM renormalization scheme and Wis 



In this section we review the basic characteristics of the RI/MOM renormalization scheme, im- 
plemented with the lattice regularization Q. After giving the essential definitions, we discuss the 
compatibility of the RI/MOM scheme with (lattice) Wis at large renormalization scales fi, where 
the presence of the Goldstone pole becomes negligible. 

We start with basic definitions. Given the quark propagator S(x± — X2',m) = {ip{xi)vp{x2)) and 
its Fourier transform S(p;m) = J dx exp(— ipx)S(x, m) (here m is the quark mass), we define two 
"projections" (i.e. traces in spin and color space) as 



T s (p; m) 



T m (p;m) 



48 



-Tr 



iTr 
12 



dS 1 {p\m) 

dS^ 1 (p; m) 
dm 



(1) 



We also consider bilinear quark operators of the form Or(x) = tpi(x)Tip2(x) where ipf{x) is the 
quark field and V a generic Dirac matrix. For definitiveness we work with two different flavors 
/ = 1,2 with corresponding masses m\ and m 2 . Specific non-singlet bilinear operators will be 
denoted as S(x), P(x) (scalar and pseudoscalar densities) and V^(x), A^(x) (vector and axial 
currents). Given the insertion of the operator Or(x = 0) in the 2-point fermionic Green's function 



Grip) = J dxidx 2 exp[-ip(xi - x 2 )} (ip(x l )OriO)^ix 2 )) (2) 
and the corresponding amputated correlation function Ar(p) = 5 _1 ip)Grip)S" 1 ip) , the projected 
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amputated Green's function Trip) 1S defined as 



r r (p) 



1 

12 



Tr pP r A r (p)] 



(3) 



where the trace is over spin and color indices and Pr is the Dirac matrix which renders the tree-level 
value of Trip) equal to unity (i.e. it projects out the nominal Dirac structure of the Green function 
Ar(p)): 



Ts 



Pp = 75 



4 7 " 



4757m 



(4) 



Everything defined so far is a bare quantity, assumed to be regularized on the lattice | We opt 
for the lattice regularization scheme with Wilson fermions. Then the corresponding renormalized 
quantities in a given mass independent renormalization scheme are formally given by 
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Z^, 1 (a^)Z r (a / u)r r (ap; am/' 



(5) 



Note that mo/ denotes the bare quark mass of a given flavor, m/ = mo/ — mc is power subtracted 

1/2 

and logarithmically divergent and m/ is renormalized. Moreover, is the quark field renormal- 
ization and Zr the operator renormalization. The functional dependence of the above expressions 
is determined by dimensional arguments and Lorenz invariance. 

The renormalization constants in the RI/MOM scheme are defined by imposing the following 
off-shell conditions in the deep Euclidean region: the wave function renormalization is obtained 
from 

r T,ip 2 / P 2 ;m 2 /p 2 



pr = pr 




lim lim Z , 1 iaa)Ty,(aa; a 2 m 2 ) 
a->0m-<-0 V ' 
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(6) 



m 



the quark mass renormalization Z m is achieved through 



T m ip 2 /p 2 ;m 2 /p 2 



P 

" 2 

m 



P 

= 



lim lim 

a— >0 m 



L z ^i a P) Z mi a P) V mia-li] a 2 m 2 



(7) 



2 This means that the integrals of eq. (^) are really sums (a 8 x ) which run over all lattice sites, labelled by 
xi, X2, etc. Also note that the /x-dependence of renormalized correlation functions is sometimes suppressed; i.e. we 
use f(p;m) instead of V(p 2 /fj, 2 ; m 2 //x 2 ). 
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and for the bilinear operator Z? is obtained by imposing the renormalization condition 



t v {p 2 /^-rh 2 /p 2 ) 



p 2 = ^2 = I™ hm o Z i/) 1 (a/i)Zr(a/i)rr(a^;a 2 m 2 ) = 1 (8) 
rh 2 = 



For any fixed scale this procedure removes the UV divergences from all Green functions and 
thus renormalization is achieved. However, renormalization conditions must also be chosen so that 
the resulting renormalized operators transform "correctly" under the chiral group; i.e. they should 
belong to the same chiral representation as the nominal bare operators. This additional require- 
ment, which ensures that chiral symmetry survives renormalization, is only true at large enough 
scales [i. At small scales symmetry violating effects due to spontaneous chiral symmetry break- 
ing in QCD should appear as AQC-D-dependent form factors in the renormalized Green functions 
^r{p 2 / /J?', rh 2 jp 2 \ Aqcd/p 2 )- Moreover, in practical simulations, usually performed at (light) non- 
zero quark mass, explicit chiral symmetry breaking form factors, proportional to m(/x) will also be 
present. Both types of form factors become negligible if the scale [i is adequately large. Therefore, 
the requirement Aqcd *C has to be satisfied in principle, so that quantities renormalized in the 
RI/MOM scheme respect chiral symmetry. Moreover, since one is usually working at fixed UV 
cutoff in simulations, the requirement [i <C 0(a~ 1 ) must be satisfied in order to control discretiza- 
tion errors. Thus, the reliability of the RI/MOM scheme on the lattice is ensured provided one is 
working in a renormalization window Aqcd <C [i <C 0(a^ 1 ). 

The above general statements have been shown to be true in ref. Q (and generalized for 
some cases of additive renormalization of dimension-six four-fermion operators in refs. |4|). In 
particular it has been shown that the RI/MOM scheme is always compatible with vector Wis, 
whereas it is only compatible with the axial Wis at large scales \i. The key observation is that the 
vector Wis 

f s (p;rh) = T m {p;m) 

f v (p;rh) = f E (p;m) (9) 

at momentum p 2 = fi 2 , are automatically satisfied by the renormalized quantities Ts, r m , Ty and 
Tg determined in the RI/MOM scheme by eqs. (^), (R) and (||). Conversely, if we use the RI/MOM 
scheme to fix, say, Ts and T m (i.e. and Z m ), then Zy and Z$ (through the identity Z$ = Z^ 1 ), 
as fixed by the Wis (^), are compatible to those obtained from the RI/MOM condition (||). 

The compatibility of axial Wis to RI/MOM is more intricate; the proof would proceed in 
exactly the same way as in the vector case, if it were not for the extra terms on the l.h.s. of Wis 



fp(p;m) + ™^^ = t m ( P ;fh) 
dm 



f a (pi m = 0) + lim t^Tt 

<?p^o 48 



757 8 " K a(p + l/ 2 ' P ~ 9/2; m = 0) 



f E O;m = 0) (10) 



which is the axial analogue to the vector Wis (^). Note that the second term of the l.h.s. of 



Wis (10) does not vanish in the chiral limit, due to the presence of a Goldstone pole. As shown in 
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ref. (U these extra terms are negligible in the deep Euclidean region Aqcd P- In this limit, the 
Wis and the RI/MOM conditions are compatible, just like the vector case. This limit corresponds 
to the lower bound of the renormalization window requirement, inherent in the RI/MOM scheme. 
Clearly, all statements made thus far are true up to discretization effects, which are present both in 
the RI/MOM and the WI implementations on the lattice. We will address this problem at a later 
stage of this work. What we need to consider at present, is that RI/MOM results are reliable in 
the renormalization window Aqcd <C p <C 0(a~ 1 ), whereas WI results are reliable in the window 

The main aim of the present work is to establish whether there is significant Goldstone pole 
contamination of the lattice RI/MOM results obtained thus far. In particular, the RI/MOM 
determination of Zp/Z$ at current (3 values and momenta does not display a clear plateau; see 
refs. J|]-@. The problem may be either due to the Goldstone pole problem or, as claimed in ref. fl2|l , 
due to discretization effects^]. Since the WI determination of the ratio Zp/Zs does not suffer from 
Goldstone pole contamination, comparison of the WI and the RI/MOM results for this quantity 



should offer a way of estimating this effect. For example, the WI result from eqs. (||) and (10) is 

Z P T P 



z s 1 + ^^£p 

Tp dm 

whereas the RI/MOM result from eq. (@) is 

¥,-¥, ^ 

The denominator on the r.h.s. of the WI (|ll|) is missed by the RI/MOM eq. (^) and this results 
to a Goldstone pole contamination of the latter. This particular WI determination, however, is 
difficult to implement in practice (e.g. realization of derivatives w.r.t. the quark mass in lattice 
numerical simulations etc.). In the next section we discuss several Wis, which are equivalent (up 
to discretization effects) but more useful in practice. 



3 The ratio Z P /Z S from Wis 

We will now review three methods, based on Wis, for the determination of the scale independent 
ratio Zp/Zg- A first method consists in computing Zp/Zg as the ratio of the PCAC current quark 
mass to the bare (subtracted) quark mass ||. In standard notation (see refs. [jOf , flltf ) the current 
quark mass for Wilson fermions is obtained from the standard PCAC relation (an axial WI, valid 
Vx ^ 0): 

9 r y dp fd 3 x{A (x)P(0)) 

2 [m - m(m )] = Z A ; ,3 x(P(x)P(0)) ( 13 ) 



3 Moreover, the 1-loop Perturbation Theory PT (and Boosted PT -BPT) result for Zp is in stark disagreement 
with the RI/MOM estimate. This comparison, however, does not reveal much, as the BPT result also suffers from 
0(<?o) errors. 
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whereas the mass power subtraction mc can be obtained by linearly extrapolating either the square 
of the pion mass m\ or [mo — fS(mo)] to their vanishing value. Then vector and axial Wis (see 
refs. [17, pi], [l6|]) determine the quark mass renormalization 

m = Zg 1 [mo — mc] = Zp 1 [mo — m] (14) 
and thus by computing [mo — m] at several mo's the ratio Zp/Z$ is obtained as the slope of 

_ Z P 

to - m = — [m - m c ] ■ (15) 
Zs 

This determination does not depend on the critical quark mass mc- It does, however, depend on 
the determination of Za- Since the methods relies on the WI (^) on hadronic states, we will label 
it Wlh. 

A second method consists in combining eq. (15) with the axial WI 

(m i + m 2 - 2m)I> (ap; ami, am 2 ) = -^Tv S~ l (ap; ami) + -^Jk S~ l (ap; am 2 ) (16) 

(considered in the mass-degenerate case), in order to obtain 

Zp Tr S~ l (ap; am) 



Zs 12 [mo — mc] Tp(ap; am, am) 



(17) 



This determination requires the projected amputated correlation functions and the quark propa- 
gators in momentum space. We anticipate at this stage, that this method fails in practice. This 
is due to the large discretization errors of the quark propagator (a detailed investigation of these 
effects, which remain large even when Clover improvement is implemented beyond tree-level, has 
been performed in ref. |12[| ). 

A third method consists in combining eq. (|l^) for the degenerate mass case with the vector WI 
(m 2 - m i)T s (ap;ami,am 2 ) = - -^Tv S^ 1 (ap; am{) + -^Tv S~ l (ap; am 2 ) (18) 



(19) 



so as to eliminate the quark propagators. Using also eq. (|15|) we obtain 

Zp (mi — m 2 ) F s (ap; ami , am 2 ) 



Zs miTp (ap; ami, ami) — m 2 Tp (ap; am 2 , am 2 ) 



The advantage of this determination is that the large 0(ap) discretization errors, present in the 
quark propagator (cf. ref. JT^[), are eliminated. Since this determination depends on correlation 
functions of external quark states, we will label it Wlq. Note that this method is characterized by 
non-degenerate quark masses. 

Finally, we mention that a WI with hadronic states has also been used for the computation of 
the ratio Zp/Zs; see refs. [16, 
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4 Results 



Our results are based on earlier simulations; all technical details can be found in refs. (2|, [l4|]. Here 
we just mention that the dataset we analyze has been obtained with the Wilson (unimproved) 
action, at couplings (3 = 6.2 and 6.4 and at several hopping parameters in the range of the strange 
quark mass. All results shown here are extrapolated to the chiral limit. The RI/MOM results are 
those of ref. B (which agree with the ones of ref. pi]) whereas the Wlh and Wlq results are new. 
We mostly show the Wilson /3 = 6.2 case 0. 

In figure [l] we show the data and linear fit characteristic of the method Wlh. The slope of the 
fitting line is Zp /[ZsZj\\. In order to extract Zp/Zg from the slope, we use the RI/MOM estimate 
for Za, as given in ref. 0. This is in accordance (with larger errors) to the RI/MOM result of 
ref. §. 



Figure 1: The value of 2[itiq — m\/ Za, obtained from eq. (|ll|), as a function of the Wilson hopping 
parameter 1/K at (3 = 6.2. The slope of the linear fit is the Wlh result for Zp /[ZsZa], according 
to eq. ( |i~5|) . The Kc value implicit in the abscissa axis is determined from the vanishing of m 2 
(linear in 1/K). 




In Fig. ^(a) we show our f3 = 6.2 results obtained from WI (|l7|), as a function of the (lattice) 
scale 

a 2 JI 2 =J2 s[n2 ( a l l v) ( 20 ) 

V 

4 Note that our (3 = 6.0 dataset has also been analyzed, with qualitatively similar results. 
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Due to enormous systematic errors, our results range over a large interval and no plateau is seen. 
The tree- level value of eq. (17) is given by 



Zp _ x , 2E,sin 2 (a Pl/ /2) 
Zs am 

where the second term on the r.h.s. is the tree-level 0{a) effect arising from the propagator sum 
in the numerator of eq. (|l7|). In Fig. ^(b) we have corrected our results by this factor, in the spirit 
of KLM tree- level improvement (cf. ref. [0]). Although this factor accounts for a big part of the 
discretization errors, the remaining effects are still large and render the method unreliable. 



Figure 2: Zp/Z$ results at j3 = 6.2. (a) Results obtained from WI fll7|); (b) Results corrected by 
the KLM factor of eq. (0). 
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In Fig. H we compare the j3 = 6.2 results for the RI/MOM and Wlq determinations, as a 
function of the (lattice) scale a 2 ~p 2 . It is clear that at small scales there is a large discrepancy 
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due to the sizeable Goldstone pole contamination of the RI/MOM method, which decreases with 
increasing scale. At large scales it is impossible to discern whether the discrepancy is due to 
a remnant Goldstone pole contamination or the discretization effects. Moreover, even at larger 
scales, the renormalization window of the RI/MOM method, compared to the plateau displayed 
by the Wlq results, is rather poor. In ref. Q, the central values of the RI/MOM results were 
taken at (ap) 2 ~ 0.8 and the systematic errors estimated as the spread of the values in the region 
0.8 < (ap) 2 < 1.5. This choice takes into account the instability of the renormalization window, 
yielding a RI/MOM result with large errors. On the contrary, the WI data displays a good plateau 
already at rather small scales. Our choice (ap) 2 £ [0.2, 0.8] for the Wlq plateau gives very accurate 
results. 



Figure 3: Zp/Z$ results at (3 = 6.2. (a) Wlq and RI/MOM determinations as a function of the 
renormalization scale, (b) Enlargement of (a), highlighting the p values where scale independence 
is expected. 
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In order to estimate the discretization errors, we collect in Fig. ||, the various determinations 
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of Zp/Zs at = 6.2 and 6.4. The most accurate estimates are those based on the two Wis (Wlh 
and Wlq). They are not, however fully compatible, due to discretization effects (recall that WI 
results do not suffer from the Goldstone pole problem). This indicates that these effects introduce 
a systematic error of ~ 15% at [3 = 6.2, which decreases to ~ 10% at (3 = 6.4. The RI/MOM 
determination, which is far less accurate due to the instability of the renormalization window, 
is nevertheless in the same range of values. Thus, it is clear that any remnant Goldstone pole 
contamination in the RI/MOM result is not superior in magnitude to the uncertainties arising 
from discretization effects. Also, the fact that the spread of the various determinations of Zp/Zs 
decrease with increasing (3, indicates that the discrepancies are due to (decreasing) finite cutoff 
effects. 



Figure 4: Zp/Zs results (a) at (3 = 6.2; (b) at (3 = 6.4. Comparison of non-perturbative 
(Wlq, Wlh, RI/MOM) and perturbative (BPT, PT) determinations. 
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Our new WI results demonstrate that the discrepancy of 10 — 15% for Zp/Zs carries over to 
Zp, computed (non-perturbatively) from 
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(22) 



In practice we use the Wlq determination for the ratio Zp/Zs in the above equation; i.e. combining 
eq. ( |j~9|) with the RI/MOM renormalization condition (|8|) for the scalar density, we have 



Zp = 



JTll — 771-2 



miYp (afx; ami, ami) — m^ p (a/i; ami, am?) 



(23) 
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Method 


Z P 


Z s 


Zp/Z s 




PT 


0.78 


0.86 


0.91 


6.2 


BPT 


0.65 


0.77 


0.84 




RI/MOM g 


0.50(5) 


0.72(1) 


0.69(7) 




Wlq [This work] 


0.57(1) 




0.79(2) 




Wlh [This work] 


0.50(2) 




0.68(2) 




PT 


0.79 


0.86 


0.91 


6.4 


BPT 


0.67 


0.79 


0.85 




RI/MOM || 


0.57(4) 


0.74(1) 


0.77(5) 




Wlq [This work] 


0.60(1) 




0.80(2) 




Wlh [This work] 


0.54(2) 




0.73(2) 



Table 1: Perturbative (PT and BPT) and non-perturbative (RI/MOM, Wlq, Wlh) estimates of 
renormalization constants. The RI/MOM values are at p?a 2 = 0.8. The Wlq and Wlh results are 



obtained from eq. (2i). 



The wave- function renormalization constant Z<j, has been computed in ref. [§]. In Table [T] we collect 
our results for Zp, Zs and their ratio. Consequently, at these [3 values, determinations of the quark 
mass (and chiral condensate) which differ in the choice of non-perturbative renormalization (i.e. 
Wlq, Wlh, RI/MOM), would also display the same variations These differences should vanish 
upon extrapolation to the continuum limit. 



It is interesting to note that, upon expressing the first of Wis (18) in terms of renormalized 
quantities and using the RI/MOM renormalization condition for Tg, the relation 



[mi(/z) - m 2 (/u)] RI/MOM 



lTr<S- 1 (//;mi)-iTr5- 1 (M; 



mi, 



RI/MOM 



(24) 



can readily be obtained. This gives a non-perturbative quark-mass determination which does not 
require knowledge of Zg (or Zp); cf. ref. [O]. Moreover, it does not suffer from large 0(ap) 
discretization effects (they cancel in the difference of quark propagators at any quark mass values), 
which would otherwise have to be isolated by chiral extrapolation. This determiantion does however 
require the computation of the wave function renormalization Z^ in the RI/MOM scheme. Since 
this is a mass-independent renormalization scheme, the above can only be implemented in the limit 
m 2 /^ 2 < 1. 
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